Bessel beams of two-level atoms driven by a linearly polarized laser field 
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We study Bessel beams of two-level atoms that are driven by a linearly polarized laser field. 
Starting from the Schrodinger equation, we determine the states of two-level atoms in a plane-wave 
field respecting propagation directions both of the atom and the field. For such laser-driven two-level 
atoms, we construct Bessel beams beyond the typical paraxial approximation. We show that the 
probability density of these atomic beams obtains a non-trivial, Bessel-squared-type behavior and 
can be tuned under the special choice of the atom and laser parameters, such as the nuclear charge, 
atom velocity, laser frequency, and propagation geometry of the atom and laser beams. Moreover, 
we spatially and temporally characterize the beam of hydrogen and selected (neutral) alkali-metal 
atoms that carry non-zero orbital angular momentum (OAM). The proposed spatiotemporal Bessel 
states (i) are able to describe, in principle, twisted states of any two-level system which is driven 
by the radiation field and (ii) have potential applications in atomic, nuclear processes and quantum 
communication. 

PACS numbers: 03.65.-w, 03.75.Be, 32.80.Qk, 41.85.Ew, 42.50.-p 



I. INTRODUCTION 

The discovery of non-diffracting light fields by Durnin 
et al. [l| triggered an interest in studying Bessel beams of 
photons @,|^. An intriguing property of these beams is 
that their intensity does not significantly change its shape 
over a propagation distance of a few centimeters [H, HJ ■ 
The Bessel beams can also have several unusual features. 
For instance, they can carry a non-zero OAM and, thus, 
make their wavefronts rotate around the propagation axis 
while the Poynting vector draws a corkscrew, so called 
twisted photons [f|. The simplest object demonstrated to 
have such characteristics is the Laguerre-Gaussian beam 
constructed in the seminal paper by Allen et al. @]. 
Twisted photons, moreover, have lead to recognizable 
advances in optical tweezers [H, Q, atom trapping and 
guiding [l(|, transfer of OAM to a system of atoms fllj ]. 
influence of OAM on beam shifts [l2[ and in other diverse 
applications (l3j ). 

In analogy with twisted photons, Bliokh et al. pro- 
posed a method to construct non-plane-wave solutions 
of Schrodinger [l4[ and Dirac [l5| equations for a free 
electron that exhibits non-zero OAM, called twisted elec- 
trons. These vortex beams of electrons have been pro- 
duced experimentally for an energy of ~ 200 — 300 keV 
only recently corresponding to the non-relativistic regime 
[l6| . The electron vortex beams have found different ap- 
plications as well, such as the magnetic mapping with 
atomic resolution in an electron microscope with twisted 
electrons [l7j and improvement of electron microscopy of 
magnetic and biological specimens [lH ■ 
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Although the photon and electron vortex beams have 
been quite extensively investigated, so far the first contri- 
bution in exploring the atomic Bessel beams, to the best 
of our knowledge, has been illustrated only in In 
an effort to extend the study of photon and electron vor- 
tex beams to other type of beams, we here demonstrate 
the construction of Bessel beams of light two-level atoms 
that are resonantly driven by the plane monochromatic 
laser beam. 

In this work, we investigate Bessel beams of two-level 
atoms that are (resonantly) driven by a linearly polarized 
laser field. For these beams, we obtain solutions of the 
Schrodinger equation in terms of the laser phase by us- 
ing an (invariant) approach as known from the rclativistic 
theory [2fj]. We utilize these solutions in order to con- 
struct atomic twisted states that go beyond the typical 
paraxial approximation. Furthermore, we examine the 
spatial and temporal characteristics of the atomic Bessel 
beam profile. In particular, detailed calculations are per- 
formed especially for probability density of hydrogen and 
light alkali-metal atom beams, that propagate perpen- 
dicular to the propagation direction of the laser beam. 
In this crossed-beam scenario, we show that the profile 
of atomic beam obtains a Bessel-squared-shape that is 
modified as the nuclear charge enlarges and can be tuned 
under the proper choice of laser parameters, such as the 
(resonant) frequency and the electric field strength. In 
addition, we also exhibit a possible enhancement of the 
second maximum in the profile of the potassium beam 
should the atom evolve in the field within a relatively 
long period, though in a time less than that of atomic 
decay. 

The paper is organized as follows. In the next sec- 
tion, we define the twisted states of a freely moving 
atom and describe their main properties. In order to 
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see how the laser beam influences the Bessel-squared- 
type distribution of the atomic beam profile, in Sub- 
section IIII A| we define the classical laser field and the 
geometry of the laser and atomic beams propagation, 
namely the collinear- and crossed-beam scenarios. For 
the interaction of a two-level atom with a linearly po- 
larized monochromatic plane-wave field, equations are 
then derived for the atomic probability amplitudes within 
the center-of-mass frame of the atom. Exact analyt- 
ical solutions to these equations are found within the 
cikonal (EA), the long- wave (LWA) and the rotating- 
wave (RWA) approximations. These solutions are then 
utilized in Subsection IIII Bl in order to construct the 
twisted states of laser-driven two-level atoms in both, 
collinear- and crossed-beam scenarios. In Section llVl we 
analyze the probability density of these (Bessel) beams 
of hydrogen and selected alkali-metal atoms, such as Li, 
Na and K, when (resonantly) driven on the Is <-> 2p, 
2s <-> 2p, 3s o 3p and 4s «-> 4p atomic transitions, re- 
spectively. In particular, we display and discuss the spa- 
tiotemporal characteristics of these beams in the crossed- 
beam scenario. Finally, a few conclusions and proposals 
are drawn in Section [V] and a detailed derivation of solu- 
tion to the Schrodingcr equation (with respecting propa- 
gation directions both of the atomic and the laser beams) 
is illustrated in Appendix |A"1 



II. TWISTED STATES OF FREE ATOMS 

A Bessel beam of any quantum particle is defined as 
a (twisted) state with its well defined energy £, longi- 
tudinal momentum pii, absolute value of the transverse 
momentum p± as well as the pro jection h£ of the OAM 
on the propagation axis [(| [l5l.l2l| . In accordance to this 
definition, the spectrum of such states can be represented 
in the form 



i>i (p) = <5(£ -£o)S(p±. ~Pxo)S(p\\ -p\\o) 



a u<t> 
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that contains a vortex phase dependency e ^, and where 
I is an integer. Such spectrum means that the momen- 
tum of the ato m is distrib uted over some cone with slant 
length po = yjp^ + p\o = const and fixed polar an- 
gle 6q with regard to the z-axis [cf. Figure [TJ the blue 
sketch]. The 'size' of this cone is defined such as the 
longitudinal and transverse components of the beam mo- 
mentum are equal to p|| = Po cos#o and p±o = po sin So, 
respectively. Owing to the conical symmetry of the mo- 
mentum distribution, we shall use cylindrical coordinates 
p = (pA.,<fi,p\\) — (psin#, (j),p cos 9) in momentum space 
and construct a Bessel-type solution of the Schrodinger 
(free wave) equation 

%l)t (r,t) = J $ t (p) *p PW {r,t)p ± dp ± dq>d Pll (2) 

as a superposition of plane waves 

ifi pw {r,t) = e^ p - r - £ V (3) 




FIG. 1. (Color online). Collinear-beam (a) and crossed-beam 
(b) scenarios for the interposition of linearly polarized laser 
and atomic Bessel beams. Red sketches indicate directions of 
the field polarization and the wave vector. Blue sketches indi- 
cate the momentum distribution of the atomic Bessel beam. 
A twisted two-level atom driven by the laser field is a superpo- 
sition of the orthonormalized waves (|18[) (a) or (|20[) (b) with 
a fixed conical momentum spread po = const, polar angle 9o 
and the azimuthal phase factor e 11 ^ . 



over the spectrum |T]) of such a cone. 

Integration ^ can be readily performed if we take into 
account the cylindrical symmetry of the atomic beam 
propagation and use the cylindrical coordinates in the 
position space as well, r = (r,(p,z). These coordi- 
nates, along with the cylindrical coordinates in momen- 
tum space, enable one to re-write the scalar product as 
p ■ r = p±r cos (0 — ip) +p\\Z in the plane wave ([3]) and, 
therefore, to exploit the integral representation 



2jt 



(4) 



of the Bessel function [22J. Direct integration simply 
leads to the final form of the twisted state of a free atom 



ifj e (r, t) = e^° z - £ot) e Ulp Ji (£) , 



(5) 



where £ = p^r/h is a dimensionless transverse co- 
ordinate which characterizes the width of the beam. 
As seen from Eq. ([5]), the state ipi(r,t) represents an 
atomic beam that propagates freely along the z-direction, 
gtpnoa/ft and nas tne p ro fii e 



(6) 



a Besscl-squared-shape in the radial dimension (as also 
illustrated in Figure [5]). The vortex phase factor e ltv 
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FIG. 2. (Color online). Distributions of probability density 
(in arbitrary units) over the dimensionless transverse coordi- 
nate for free (gray solid curve) and laser-driven (blue dashed 
curve) beams of twisted atoms with OAM M — 2ft. For the 
crossed-beam scenario, the probability density is depicted 
at t = when the field is switched off. The distribution is 
shown by the variation of colors from black to white within 
the 'sunset' scale (left panel), where black and white corre- 
spond to the minimum and maximum values of the probability 
density. 



ensures that the twisted state (|5|) is an cigenstate of the 
z-component of OAM operator £ z = —ihd/dip, and, thus, 
is responsible for non-zero OAM M on the propagation 
axis. In addition, one should stress that the states (O 
are orthogonal and can be normalized if the integration 
is carried out, for instance, over a large, but finite cylin- 
drical volume. For our further analysis, however, the 
normalization is not of crucial interest. 

We have shown a simple, but very significant procedure 
how one can create twisted states of a free particle (e.g. 
atom) . In our later study, we are interested in how these 
twisted states are affected by an external influence, such 
as a driving laser field. For this purpose, we will initially 
prepare a two-level atom in the upper state and then 
switch on the laser field. Once this atom is coherently 
driven by the field, we will apply the same approach in 
order to construct Besscl beams of such atoms and to 
investigate their spatiotemporal characteristics. 



III. TWISTED STATES OF LASER-DRIVEN 
ATOMS 

In the previous section, we have shown a simple way 
of construction of free atomic Bessel beams. In this sec- 
tion, we are interested in the dynamics of these beams 
when they are (coherently) driven by the field of a plane 
monochromatic laser wave. To this end, in subsection 
IIII Al we examine the coherent coupling of a two-level 
atom with the linearly polarized radiation field, and then, 
in Subsection IIII B[ we construct the twisted states of 
this laser-driven atom by employing the (Bessel) spec- 
trum (fTj). 



A. Semi-classical coupling of two-level atoms to a 
laser field 

The coherent interaction of a two-level atom with the 
radiation field has been explored since the early days 
of quantum optics [23hj25j . Nowadays, this atomic co- 
herence is known as an effective tool for achieving con- 
trol about (atomic and molecular) samples [2fjj and it 
can be utilized to examine various effects with atomic 
beams, such as generation of entanglement f27j . investi- 
gation of the optical force [28| or exhibition of vortices 
[29( in beams. In order to examine how this coherent 
control can be exploited also for tuning and guiding the 
atomic Bessel beams, in the following, we shall describe 
the coherent coupling of an atom to the laser field in the 
meantime respecting the propagation directions of both 
the atom and laser beams. 

To describe the interaction of a two-level atom with a 
classical radiation field and to start our derivations, let us 
first characterize the field and the interaction (Hamilto- 
nian) between the atom and this field. We hereby assume 
that the atom has the (rest) mass m, a dipole moment 
d, and that it moves with constant momentum p along 
some given direction. 

1. Characterization of the classical field 

We consider the classical electric field of laser beam as 
a (monochromatic) plane wave 

E = £e i(k-r- U t) 5 ( ? ) 

with constant electric-field amplitude e, frequency ui and 
wave vector k = kn, which satisfies Maxwell's equa- 
tions, and where the unit vector n defines the propa- 
gation direction of the wave [3(|. While, in general, E 
and £ are both complex- valued vectors in Eq. ([7]), the 
physically relevant electric field is given by the real part, 
$t(E), and is for transverse waves always perpendicular 
to its propagation, i.e. n ■ e = 0. For our further dis- 
cussion, moreover, it is useful to introduce a set of real 
and mutually orthogonal unit vectors (ei , e% ,n), and to 
re-write the field amplitude in terms of these vectors as 

£ = ei s + e 2 e' , (8) 

with the two (complex) constants Eq and e' . 

Most generally, the laser field amplitude e in Eq. ([7]) 
describes an elliptically polarized plane wave; as a spe- 
cial case, this definition includes linearly polarized waves 
if the complex constants £o and e' fulfill proper relations. 
We choose e' Q = which corresponds to a linearly polar- 
ized field in e± direction 

E(r,t)^e 1 e Q e l{k - r - ult) 

with the (real) amplitude £o- Therefore, for this wave 
the relevant electric field is given by 

E = ciEq cosC = ecos(, (9) 
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where C = k ■ r — uit is the phase of the plane wave 
independent of its particular polarization properties. We 
here should stress that the field is evaluated at the center- 
of-mass position r = (x,y,z) of the atom. 

The field and the atom propagation directions, in prin- 
ciple, can be chosen arbitrarily which will form the ge- 
ometry of the "atom + laser" system. To describe the 
linearly polarized laser field we specify this geometry al- 
ready in this subsection. We examine two scenarios that 
depend on the interposition of the propagation direc- 
tions of laser and atomic beams. These are the so-called 
collinear- and crossed-beam scenarios for which the laser 
and the atomic beam propagation directions correspond- 
ingly are parallel or perpendicular to each other [cf. Fig- 
ure!]]. As mentioned above, the z-axis is chosen along the 
propagation direction of the atomic beam. In collinear- 
beam scenario, moreover, we chose the x- and z-axes di- 
rected along the polarization t\ and laser beam propa- 
gation, respectively. In contrast to this, in crosscd-beam 
scenario, we restrict ourselves with the polarization along 
z-axis and consider propagation in x — y-plane, where the 
x-axis is declined from the laser propagation direction 
under (jj^ angle. 

With the distinction of these collinear- and crossed- 
beam scenarios the laser field © and the phase £ obtain 
the form 

E {ll) = (ecosC (ll) ,0,0), C (ll) = k\\z-ut, (10) 
E (±) = (0,0,ecosC U) ) , C (±) = kj_rcos(<j> L - <p) - ut, (11) 

respectively, with a transverse coordinate r = \/ x 2 + y 2 
of the atomic center-of-mass. In the following, we define 
the Hamiltonian of our system and determine the atomic 
states that are driven by the laser field (fit)]) and (fTTj) . 



2. Solution of the Schrddinger equation in the 
center-of-mass frame of the atom 

In the semi-classical theory, the quantum dynamics of 
the atom is driven by the external field but this "motion" 
does not re-act back upon the field as this would require a 
quantization of the electromagnetic field (23|. Therefore, 
the Hamiltonian of an atom in the (classical) field takes 
the form 

H = 1- -ffatom + -Hint, (12) 

2m 

"2 

where w— denotes the kinetic energy (operator) as asso- 
ciated with the center-of-mass motion and 

H atom = E a \a){a\ + E b \b){b\ (13) 

refers to the internal motion of the atom, whilst the in- 
teraction of the atom and the field is taken (as usual) in 
minimal coupling 



in the LWA [25[. The two state vectors \a) and \b) in 
expression (| 13[) denote the upper and lower states of the 
atom and are supposed to be eigenstates of the atomic 
Hamiltonian 

^atom I a) = E a | a) 
ffatom \b) = E b \b) , 

and where E a and E b are the energies of upper and lower 
states, respectively. 

Making use of the completeness \a) (a\ + \b) (b\ = 1 for 
a two-level system, the interaction Hamiltonian of the 
atom with a linearly-polarized wave can be written in 
the form 

= -m^e^ \b) (a\ +e- 4 **» |o) (&|)cosC ( "\ (14) 

= -ttlH.(e* d - \b) (a\ + e-«* \a) (6|)cosC (i) (15) 

for collinear- and crossed-beam scenarios, respec- 
tively. In accordance to these scenarios, here 0,r x = 
|e (b \x\ a) | £q/H and fi^. = |e (b \z\ a)\ So/h denote the 
Rabi frequencies that describe the coupling strength of 
the atom either with x- (fTOj) or z-polarizcd (fTTj) field, 
respectively. With this distinction for the Rabi frequen- 
cies, the exponentials 4>d T and 4>d, denote the phases of 
the dipole matrix elements: e (b \u\ a) = \e (b \u\ a) |e^ d " 
where u = x, z. Here both x and z refer to so-called inter- 
nal variables of the atom and are the electron's relative 
coordinates with regard to the nucleus [25| . The Rabi fre- 
quencies, moreover, arise naturally as non-diagonal ma- 
trix elements of the atom dipole moment if wc express 
the interaction Hamiltonian in the {|o) , |6)} basis. This, 
in turn, makes the interaction Hamiltonian dependent on 
atomic states |a), \b) and the center-of-mass coordinate 
of the atom [3lj] . 

For such an effective factorization of the interaction 
Hamiltonian, and as described in many texts before 
[23l . [24| , the two- level atom undergoes Rabi oscillations 
between its lower and upper states with frequency Qr^ 
or fln z , quite analogue to a spin- 1/2 system in an oscil- 
lating magnetic field [32|. In contrast to the standard 
derivation (see, for example, Ref. [23[), however, the in- 
teraction Hamiltonian (Ti~4")) - (|15p now also depends on the 
phase and ( {±) of the radiation field to account for 
the time- and space-dependency of the atom-laser inter- 
action. In addition, some time ago an approximate tech- 
nique, i.e. expansion of the exponent (J7J in the inter- 
action Hamiltonian, has been applied to the helium and 
hydrogen atoms in order to examine their coupling to an 
intense laser field [HI]. 

To explore the time evolution of the atom, let us search 
for solutions of the (time-dependent) Schrbdinger equa- 
tion 

Hib = ih^$- . (16) 
dt x ' 

If wc utilize the (effective) factorization of the interaction 
Hamiltonian we can use the ansatz 



Hi n t — d ■ E 



1> (r, t) = e 4(f— f *) (V> a (C) \a) + MO |6» , (17) 
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as a solution of the Schrodingcr equation. Here £ = 
( {±> is the laser phase for collinear- or crossed-beam sce- 
narios, respectively. The constant (non-relativistic) mo- 
mentum p = mv of the center-of-mass of atom gives 
hereby rise to the space-dependent "translation factor" 
en pr to account for its overall motion with energy £ 
within the given coordinates. A similar ansatz has been 
applied in [34| in order to investigate Stark splitting of a 
three-level atom. 

We bring the detailed derivation of solution of the 
Schrodingcr equation (|16[) in Appendix|X]where the phys- 
ical assumptions we consider are the following. We pre- 
pare the atom initially in the upper state [cf. Eq. (|A16I) ] 
and then employ the typical, eikonal and rotating wave 
approximations, for which the atom rest energy is much 
larger than the photon energy (EA) that, in turn, is res- 
onant to the atomic transition energy (RWA, see also 
Eq. (|A17[> ). For light two- level atoms,both the EA and 
RWA are valid with high accuracy J23(. Under these 
physically relevant conditions, we obtain analytically ex- 
act solution 

^(11) ( Pjt ) = e i( P -r-£t) ( e -ia<MC™ COS-^- \a) 

+ te^^e ip c sin ^ — | £,) j (18) 

that represents the state of the laser-driven two-level 
atom for collinear-beam scenario with some reduced 
quantities [cf. Eq. (|A3[) ] 

{ a (H> , ffdi) , Q0i>} = {Ea ,E b ,n B } ^ 
h [kvn — LOj 

Whereas, in crossed-beam scenario, the atomic state is 
given by 

^ (r,t) = e ^-zt) ( e -^ [±} C^ cos 9^1 \ a ) 



with reduced quantities 

h{a (± \ f3 (± \ n^}- 



kv± cos(<Pl 



(21) 



In Eqs. (fI5 ]) -(|2T ]) . p 2 = pi + p\ is the squared to- 
tal momentum of the atom (with its longitudinal p\\ 
and transverse p± components and the atomic veloci- 
ties U|| = p\\/m and v± = p±/m, respectively). The 
orthonormalized solutions (fT8| and (|20j) represent a su- 
perposition of the upper \a) and lower \b) states with £ (ll> - 
or C (J - ) -dependent coefficients. The 'translation' factor 
e~fi(P' r ~ £t \ as mentioned above, describes the motion of 
the atom as a whole with momentum vector p along some 
(chosen) direction and energy E. 

So far, the focus of this subsection was placed on de- 
termining solutions of the Schrodingcr equation (|16|) for 



linearly polarized field (p~0|) or (JTTJ) . In the next subsec- 
tion, we will utilize the explicit r-dependency of (exact) 
solutions (fT8"|) and ([2"U]) and exploit them in order to con- 
struct a Bessel beam of laser-driven two-level atoms. 



B. Twisted states of laser-driven two-level atoms 

In the previous subsection, we have analyzed the timc- 
and space-dependent interaction of two-level atoms with 
a linearly polarized laser field and have built the states 
(fTB"]) and (J2T)|) which describe the spatiotemporal dynam- 
ics of laser-driven atoms. In the following, we shall 
utilize these states in order to construct atomic Bessel 
beams which carry a non-zero OAM and meanwhile do 
not diffract along the propagation direction. 

The wave function of the laser-driven two-level atom 
with the projection M of the OAM on the (atomic beam) 
propagation axis can be constructed as a superposition 
of orthonormalized solutions (fTB"]) and (f2T)|) of the time- 
dependent Schrodinger equation (ll(i|) 



^l(r,t) = J ip t (p) ip(r,t)p^dp ± d(l}dp ll (22) 

over the monoenergetic cone ([I]). This step in the con- 
struction of twisted atomic beams is analogue to the use 
of twisted photons [f| and electrons fl5j | . 



1. Collinear-beam scenario 

Let us start with the collinear-beam scenario for which 
the integration (|22p can be carried out easily in a similar 
way as for the twisted state (O of a free atom. Indeed, 
by performing the same steps, i.e. substituting the state 
(fT8"| into Eq. (|2"2"|) and making use of the integral repre- 
sentation (j4]) of Bessel function, we obtain a simple form 
for the twisted state of a laser-driven two-level atom 

/ nn fin o (ll V<H' 
= e *(Piio*-£ot) e «*j,(£) ( e -^o m C m C08 l^— \ a ) 

_|_ ,. P l <Pd x p Vo s sin " 



|6) 



(23) 



Here the dimensionless transverse coordinate £ = p\_or/K 
describes the width of the beam, and the reduced quan- 
tities ceo", /3o M> and fig are taken on the cone 'surface' 
p = Po [cf. Eq. (Jin]) and Figure QJi] . 

The state (|23|) is the eigenstates of the z-component 
of OAM operator £ z due to the presence of the vortex 
phase factor e liv . Hence, apart from the free propa- 
gation along the z-direction, e lp ^° z ^ h , and the Bessel- 
dependency on the (dimensionless) transverse coordinate 
the twisted state , F^" > also carries a well-defined OAM M 
quite similar to the scalar [13, [H[ and spin-dependent 
[15j electron Bessel beams. 

It is obvious that when the field is switched off, i.e. 
fjCli) _ £01) — o, the free twisted state is recovered. 
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The presence of the field causes only Rabi flopping of 
(already) twisted two-level atoms between the upper \a) 
and lower \b) states and meantime does not influence the 
Bessel-squared-shape of the beam profile. Alluding to 
Section HVl let us note that the absence of the field's im- 
pact on the beam profile is quite expected since the laser 
phase C <M) = kuz — uit in collincar-beam scenario con- 
tributes only in the free motion of the atomic beam along 
the z-axis and does not contain the radial coordinate r in 
contrast to the laser phase in crossed-beam scenario 
where both r and the azimuthal angle ip are involved [cf. 
Eqs. QlO p -Qli p ]. It is the presence of r and <p that enables 
one to explore an intriguing dynamics of atomic twisted 
states in crossed-beam scenario as we exhibit in details 
in our forthcoming discussion. 



2. Crossed-beam scenario 

So far, we have investigated the twisted states of free 
and of laser-driven atoms in the collinear-beam sce- 
nario (|23p . These two states involve the same, station- 
ary Bessel-dependency Ji{£) in the profile of the beam. 
Whereas, the crossed-beam scenario gives rise to a time- 
dependent atomic beam profile. To investigate this we 
shall evaluate the integral ([22]) by employing the solution 
(|20p . This integration cannot be carried out by means of 
exact methods since the reduced quantities a (±) , (3 {±) and 
£l (±) depend on the </> angle [compare Eqs. (fT^) and l|2ip]. 
However, if we consider a non-relativistic regime for the 
propagation of atomic beam we can easily proceed ana- 
lytically. For this purpose, we again represent the scalar 
product in the form p ■ r = p±r cos {4> — ip) + p\\z and 
evaluate the integral (f22|) first with respect to transverse 
p± and longitudinal p|| components of the (linear) mo- 
mentum 



>te(piio*— £ot) 



x / d^e^e 



where the state vectors 



5>"-nS»>, (24) 



71 = 1 



|Bi,a> = 



B 3 a) = ± 



,i<Pd z 



\b) 



refer to the upper and lower levels of the atom. For the 
sake of brevity, moreover, we introduce the dimensionlcss 
parameters 



-^±0^/2,(25) 



which carry information about the internal structure of 
the atom and the strength of the atom-field coupling. 
These parameters are given in terms of the reduced quan- 
tities , and f2g L) and are taken at the cone surface 
P = Po- 

Now we apply our assumption about the non- 
relativistic regime for the propagation of atomic beam 



and Taylor expand A n (n = 1..4) by keeping only the 
terms up to the first power in v±q/c 



vA n ~ C n ( 1 



■ COS(0 L - 



where 



Ci, 2 = (E a T m R j2)/(fiw) , C 3 , 4 = (E b T m R j2)/(fvuj) 

are dimensionless constant energies normalized to the 
photon energy. This expansion enables one to re-write 
the wave function (|24p in a desired form 



x 

n=l 



j giXn coscf) g*3^n sin <fi 



\B n ) (26) 



which is appropriate for analytical integration. Here the 
dimensionlcss coordinates 

X n (£, <p,t)=£ cos ip + C n ( {±) — cos <f) L , 

c 

y n (£,<p,t) = £ sin <p + C n C (±) — sin L (27) 

c 

contain the (time-dependent) laser phase 
ftfc 

C ±} = £cos(0 L - ip) - wt, k = k±_ (28) 

P±o 

that is independent of the integration variable <j). 

In order to calculate the integral (|26[) we introduce a 
new system of 'cartesian' R n = (X n , y n , Z) and 'cylindri- 
cal Fin = (t_. n , $„,Z) coordinates with standard trans- 
formations 

X n = E n cos $„ , y n = E n sin $„ , Z = p± z/ ft , (29) 



where S n = yj X% + y% is the 'radial' coordinate and 
<&„ is the 'azimuthal angle' in the (X n ,y n )-pl&ne. By 
inserting these new coordinates in Eq. (|26|) and making 
use of a simple trigonometric relation, we can exploit the 
integral representation of the Bessel function 



2 71 



s iH n cosW-*„) = 2 7 ny M, "^(S„) 



and obtain the final form of the twisted state of laser- 
driven two-level atoms in crossed-beam scenario 

*w = }_ e i(rn°*-£°t) ^e iC ^ W e^J,(S„) \B n ) . (30) 

n=l 

It is important to note that we recover the free twisted 
state ([5]) if we switch off the laser field (i.e. ( {±) —> 0), 
since in this case the coordinates S n and <!>„ coincide with 
£ and <p, respectively, as (X n ,y n ) -> (p ±0 x/h,p ±0 y/h) 
[cf. Eq. (JUl)]. 

Since no restriction has been made on the longitu- 
dinal p||o and transverse components p±o of the atom 
momentum, the states (f2"3"P and (f3T))) apply generally for 
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scalar Bcsscl beams of two-level atoms beyond the typ- 
ical paraxial approximation for both, the collinear- and 
crossed-beam scenarios. As we stressed above, we are 
not interested in the normalization of these states since 
it does not provide too much insight in the overall dy- 
namics of the beam. To give a hint, however, we would 
like to point out that these beams, can be normalized, 
on one hand, by integrating them in a large, but a fi- 
nite cylindrical volume, as already mentioned above and 
also done in [2l[ for twisted photons. On the other hand, 
we could regularize the Dirac ^-function, which arises af- 
ter the integration of squared Bessel functions over the 
whole space, by rigorously re-defining it as a limit of reg- 
ular functions (e.g. Gaussian), as also mentioned in [15( 
for electron Bessel beams. 

In the last discussion of this section, we put the em- 
phasis on finding an integral of "motion" that describes 
the propagation of the laser-driven twisted atom in the 
crossed-beam scenario. A similar operator description of 
angular momentum properties of light Bessel beams has 
been done in [36|. Given that the photon energy and 
momentum are much less than the atom rest energy and 
transverse momentum, respectively, i.e. huj/(mc 2 ) <C 1 
and hk/(mv±o) < 1. wc replace (|27|) with the following 
approximate relations 

Xn ~ £ cos ip — C n v±okt cos (f>L , 

y n w £ sin ip - C n vj_ kt sin <f> L , (31) 

that are valid with high accuracy for a wide range of fre- 
quencies, from infrared to ultraviolet regions. The wave 
function (|30|) . therefore, can be re- written as 



*r 5 ~ *< B) > 



(32) 



$W = _ e *(P||o*-(&+c„J»u)t) e *»» i £(E f ,) \B n ) , (33) 



where we have replaced the laser phase ([28]) with m 
—tot. One can immediately recognize that the exponen- 
tial e lp ^° z / h describes the diffraction-free propagation of 
the beam along the z-axis. Moreover, this free propaga- 
tion occurs as a sum of four scalar Bessel modes \E'^" ) , 
each of which carries a quasi-energy £q + C n hio because 
the atom is dressed by the field. 

The notable difference between the twisted states of 
laser-driven (also of free) atoms in collinear- and crossed- 
beam scenarios is that here the Bessel function depends 
on the new coordinate S„ and, consequently, on time 
(compare Eqs. ([23]). (|5j} with Eq. ([32])). Due to this (dif- 
ferent) dependency, the state (|3"2"j) is no longer the eigen- 
state of the conventional OAM operator £ z = —ihd/d(p, 
instead each of the modes (f3"31 represents an cigenstate 
of the operator 



Ci n) = -ih 



(34) 



an OAM operator which acts in (X n ,y n , Z) configura- 
tion space, upon which our physical system depends. It 



is easy to verify that this operator has the eigenvalue M, 
i.e. = M^"'. Whereas, as mentioned above, the 

states (J3J and ([2"3"]) are the eigenstates of £ z with the same 
eigenvalue M, i.e. i z ip t = M^i and £ z ^° = M%f\ 
respectively. In addition, as one may expect the opera- 
tor (|M|) coincides with the operator £ z when the field is 
switched off, C ( z n) — > l z , (see also Eq. ([55} ) . 

To get a deeper insight let us now calculate the mean 
value of the conventional OAM operator averaged over 
the modes ^7* . For this purpose, we express £ z in terms 
of coordinates ([29]) by employing the relations (|3"Tj) 



-ih 



d_ 

dip 



-ih 



d$ n d . dE n d 



-ih 



C n v± kt 



dip d$ n dp dE. 
cos($„-0l) d 



d<P n 



C„u_L fctsin($„ - 4> L ) 



d 



(35) 



For the (quantum mechanical) mean value of £ z , when 
normalized to the overall 'energy' of the mode (1331) . wc 
obtain 



(4) 



0lRn 



hi. 



(36) 



dRn *i n 



where dJR n = S„c?S„d$„c?Z is the elementary cylindrical 
volume in configuration space which is related to the ele- 
mentary volume in position space via dr = (h/p±o) 3 dR n . 
The last two terms of the operator (|35|) do not con- 
tribute to the integral after the integration over 
As seen, the mean values of both OAM operators co- 
incide, (C z n) ) — (£ z ). This means that, apart from 
the diffraction-free propagation along z-axis, the state 
([3"2]) describes a superposition of four modes (|33| each of 
which carries a non-zero OAM with respect to the same 
axis. Moreover, the time-dependent profile of the Beam 
reveals a non-trivial dependency on the transverse coordi- 
nate £ and the azimuthal angle ip. In the next subsection, 
we examine this non-triviality in details and spatiotcm- 
porally characterize the Bessel beams of two-level hydro- 
gen and selected alkali-metal atoms that are resonantly 
driven by the laser field without the level damping. 



IV. SPATIAL AND TEMPORAL 
CHARACTERIZATION OF ATOMIC BESSEL 
BEAMS 

So far we have built the twisted states (|2"3"]l and (|3H|) . 
(f3"2")) of laser-driven two-level atoms for collinear- and 
crossed-beam scenarios. These states can be used in or- 
der to study the space- and time-dependent profile of 
atomic beams. In this section, therefore, we put the em- 
phasis on the crossed-beam scenario and investigate how 
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the radial distribution and the time evolution of the prob- 
ability density are affected by the atomic beam velocity 
and the nuclear charge Z when the laser field is resonant 
to the atomic transition energy. 

To proceed with our further discussion, we define the 
probability density of atoms in a twisted state with a 
non-zero OAM M as 



I*/ 



(37) 



such that p = 1 for beams (fTSj) and ([20]) of non- 
twisted atoms. Depending on which of the two scenar- 
ios occurs, the probability density (f3"T)) acquires different 
forms. To show this, we substitute the wave function 
([23]) into Eq. (|37|) and obtain the probability density for 
the collinear-beam scenario 



Pi 



(38) 



as a function of only the dimensionlcss transverse coor- 
dinate £ . Figure [2] shows the 'non-diffracting' distribu- 
tion of this probability density that coincides with the 
beam profile of free twisted atoms (even if the laser field 
is still switched on). This coincidence is due to the £- 
independent laser phase (| 10[) as also pointed out in Sub- 
section uttbt] 

Let us now calculate the probability density in crosscd- 
beam scenario. For this purpose, by inserting the state 
(f3"D)) into Eq. (|3T|). after straightforward derivations, we 
obtain 

pf > = et + A e , (39) 
where the (space- and time-dependent) term 
I 4 

n=l 

is responsible for the Bessel-squared-type distribution of 
the probability density, and the term 

- \ C0S (^C± + ^($4 - $ 3 )J 4 (S3) Ji (S 4 ) ,(41) 

is a small summand that can be neglected with high ac- 
curacy under properly tuned parameters of the "atom + 
laser" system, as shown below. 

In order to explore and exhibit the temporal and spa- 
tial characteristics of atomic Bessel beams let us con- 
sider, for example, two-level hydrogen, lithium, sodium 
and potassium, and assume that these atoms are driven 
on the Is O 2p, 2s f-> 2p, 3s -H> 3p and 4s •*-> Ap atomic 
transitions, respectively. For the sake of simplicity, how- 
ever, we here also suppose that no decay occurs for upper 
levels and, thus, that no damping applies in the time- 
evolution of the probability amplitudes (|A18|) - (|A19I) . For 
the 2p<-> Is transition, the laser and Rabi frequencies for 
hydrogen can be easily expressed as 

u = Sl R . = *^°, n B ,=V2fi H . (42) 
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FIG. 3. (Color online). Distribution of probability density 
(in arbitrary units) as a function of dimensionless trans- 
verse coordinate £ = p±or/h (left panel) and azimuthal angle 
ip (right panel) for i — 2. The probability densities are shown 
for Bessel beams of H (al-2) with atomic velocity 1.4 ■ 10 
cm/s and of Li (bl-2), Na (cl-2), K (dl-2) with atomic veloc- 
ity 0.7 ■ 10 6 cm/s for different propagation times 10 fs (gray 
solid curves), 15 fs (blue dashed curves) and 20 fs (red dot- 
dashed curves). 



by using the well-known properties of hydrogen-like ions 
|37J. For alkali-metal atoms, we make use of the known 
values of their spectrum [38| and dipole matrix elements 
}39l ]. Moreover, we here restrict our discussion to low- 
and medium-^ atoms since the interaction Hamiltonian 
(fI3 ]) -([T5 |l is valid only within the LWA, i.e. when the 
radiation wavelength exceeds the atomic sizes. 

After we have specified the type of a laser-driven two- 
level atom we are ready to explore both, the spatial and 
temporal characteristics of atomic Bessel beams. When 
the laser field is switched off, i.e. for t = 0, the both 
curves for collinear- and crossed-beam scenarios coincide 
[cf. Figure^]. This is quite expected since in the absence 
of the laser radiation only a free Bessel beam propagates, 
as one might also expect due to the initial conditions 
(|A16[) . Once the laser is switched on, the atomic beam 
starts evolving in the field and changing its conventional 
Bessel-squared-shape. Figure [3] displays the probability 
density profile (i) at various dimensionlcss transverse co- 
ordinates ranging from to 15 at given azimuthal angle 
tp = 7r/3 (left panel) as well as (ii) at various azimuthal 
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FIG. 4. (Color online). 'Snapshot' of atomic Bessel beam 
profiles at t — 20 fs after the laser, which drives the atom, 
is activated. The laser wave propagates along the positive 
direction of t/-axis, the polarization vector points toward the 
observer. The parameters of the "atom + laser" system and 
the variation of colors are the same as in Figures [3] and [2] 
respectively. 



angles ranging from to 2ir at given (two different val- 
ues of) transverse coordinate £ = 3 and £ = 5 (right 
panel). In particular, results are shown for the nonparax- 
ial atomic beam with transverse momentum p±o = po/5 
and OAM M = 1% as well as for the field strength e — 4 
GV/cm and laser propagation angle <pL = 7r/2 for three 
different evolution times, 10, 15 and 20 fs. Moreover, Fig- 
ureS]combines both the £- and (^-dependencies and shows 
the actual profile of Bessel beams of the same atoms when 
they propagate 20 fs in the field. 

As seen in Figure [3] the deviation of curves from each 
other increases the longer the atom propagates in the 
laser field. This deviation is caused by all four Bessel 
modes in Eq. (|40|) containing the arguments 



S„ = \ U 2 + 2£C„C< 



Vj_ 



■ cos ( 



v/£ 2 - 2£C n v±nktcos (<j> L - ip) 



(43) 



The time factor 2^C n v±okt cos ((j)L — <p\ which involves 
both the transition energy and the atom-laser coupling 
strength, can lead to an enhancement of the second maxi- 
mum of the beam profile during its evolution in the laser 
field. This enhancement depends crucially on how the 
parameters of the "atom + laser" system are tuned. For 
instance, when the beam of potassium is evolved in the 
field within 20 fs, the first two maxima become of the 
same order [cf. Figured] (dl)]. This can be clearly seen 
also (i) from Figure |3] (d2) in the region where the red 
dot-dashed curves intersect for £ = 3 and £ = 5 and (ii) 



FIG. 5. (Color online). Distribution of for I — 2 as a func- 
tion of dimensionless transverse coordinate £ = p±or/h for hy- 
drogen with different propagation times (left panel) and for 
alkali-metal atoms with an evolution time 10 fs (right panel). 



from Figure [3] especially for the azimuthal angle ip ~ ir/2 
for which the first two maxima are well separated (white 
areas for the potassium profile in the vertical direction). 

Figure [3] exhibits also another intriguing feature of 
the beam profile in the crossed-beam scenario (which, 
in fact, eventually leads to the enhancement of the sec- 
ond maximum). Since the term C n v±okt reveals different 
values for different atoms, this Z-dependency gives rise 
to field-induced 'spread' of the Besscl-squarcd-shape of 
the atomic probability density [cf. Figure [3] (bl)-(dl)]. 
This behavior is rather universal and remains the same 
for all z, i.e. along the atom propagation axis. Indeed, 
as the nuclear charge increases the white areas in Fig- 
ure 2] start to symmetrically spread along and against 
the propagation direction of the field, meanwhile keeping 
their shape constant along the z-axis. All these non- 
trivial spatiotcmporal characteristics are caused by the 
coherent interaction of atomic and laser beams and are 
quantitatively reflected both in the Z-dependent atomic 
transition energy and the atom-field interaction strength, 
i.e. the Rabi frequency. In addition, we would like to note 
that the term Ai does not contribute in the radial distri- 
bution of the probability density because it is effectively 
zero for an evolution time varying in 0, ..100 fs and for 
low- and medium- Z atoms [cf. Figure [5]. 

Asjfor example, in the case of optical [1(1 E3 and elec- 
tron [41[ Bessel beams, we have spatially and temporally 
characterized the atomic Bessel beams that are driven 
by the laser field. In contrast to these studies, radiation- 
assisted atomic Bessel beams acquire a special type of 
behavior as illustrated in Figures [3][S] due to the coherent 
coupling of an atom to the laser field. 



V. SUMMARY AND CONCLUSIONS 

The twisted states of laser-driven two-level atoms have 
been built and investigated, especially, in crossed-beam 
scenario when the laser and atomic beams are perpen- 
dicular to each other. In more detail, the interaction of a 
two-level atom with a linearly polarized electromagnetic 
field has been described by using the space- and time- 
dependent laser phase for solving the Schrodingcr equa- 
tion in a similar way as known from relativistic quan- 
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turn theory of electron. Exact analytical solution to 
the Schrodingcr equation was found within the eikonal, 
rotating-wave and long-wave approximations (to deal 
with fields nearly resonant to the two-level excitation en- 
ergy and with wavelengths larger than the atomic size). 
Our treatment enables one to construct a twisted state 
of laser-driven two-level atoms with their well defined en- 
ergy, transverse and longitudinal momentum components 
as well as the projection of the orbital angular momen- 
tum along the propagation direction. By making use of 
these states, detailed calculations have been performed 
for the distribution of the probability density of hydro- 
gen, lithium, sodium and potassium for the Is <-> 2p, 
2s -H> 2p, 3s <-> 3p and 4s -H> Ap atomic transitions, 
respectively, without (level) damping. For the crossed- 
beam scenario, we have exhibited a non-trivial, Bcsscl- 
squared-type behavior of the beam profile that applies 
for both, paraxial and nonparaxial regimes and depends 
on time. We have also shown that a possible enhance- 
ment of the second maximum of probability density may 
occur under a specific choice of laser and atom param- 
eters, such as the atom evolution time, nuclear charge, 
atomic velocity, the laser frequency and the electric field 
strength. 

Though emphasis was placed on the interaction of two- 
level atoms with external (monochromatic) fields, the 
theory in this work is applicable also for three- and 
multi-level atoms of different configurations, including 
A—, V— and S— type atoms. Moreover, this theory can 
be extended quite easily to a (general) elliptically polar- 
ized monochromatic field as well as to a standing wave 
E = (ecos(fer - ojt) + ecos(fcr + wt),0, 0). The latter 
case will enable one, for instance, to deal also with elas- 
tic interactions between the atom and the field. In the 
future, moreover, it seems desirable to us (and possible) 
to construct a unitary operator of the form 



U 



cos(r>[ ) ll) C (ll) /2) ia (io u(t _ t/) 



cos(n< ) ll) C ,(ll) /2) 

Bin(n<">c'"V2) e0 »» w( «-o , b){b] 

with the phase = fen z—uit' at some initial state in the 
collinear-beam scenario, in order to describe collision and 
scattering phenomena with atomic Bessel beams. Such 
evolution operators would be useful also for investigation 
of form factors of Bessel beams of radiation-assisted two- 
level systems, such as molecules, atoms or even nuclei 

m. 
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Appendix A: Derivation of solutions of the 
Schrodinger equation 

To determine the probability amplitudes ip a and ipb for 
the atomic states | a) and | b) , we take into account that 
the overall dynamics of the atom-laser system remains 
the same for the collinear- and crossed-beam scenarios 
(compare Eqs. (fl"4]) - (fl"5j) ), though they imply different 
coupling strengths and phases in the interaction Hamil- 
tonian. For these similar Hamiltonians, in our further 
calculations we will adopt generic notations £Ir, <f>d and 
£ in order to replace the Rabi frequencies Qr^. , Q r z , the 
dipole matrix element exponentials 4>d x , 4>d z and the laser 
phases C||, Cj-j respectively. The replacement of these 
laser phases is justified since the solution which we con- 
struct does not depend on any preferred direction of the 
laser wave vector. 

In order to describe the dynamics of laser-driven two- 
level atoms we take into account that the laser phase £ 
contains both, the time- and space-variables and, there- 
fore, enables one to express the corresponding partial 
derivatives by the total derivative with regard to this 
phase. Based on this mathematical trick we can re- 
write the Schrodinger equation as an ordinary differen- 
tial equation in a similar way as the Dirac equation has 
been examined earlier by Volkov [2(| and Skobelev [43| 
who found solutions for the (relativistic) motion of elec- 
trons and neutrons, respectively. Thus, by making use of 
this technique and substituting ansatz (|17[) into the time- 
dependent Schrodingcr equation (|16[) . we obtain the two 
coupled equations 



h 2 k< 

2m 

h 2 k' 
2m 



ip'a +ih(v ■ k — w) ip' a - E a ipa + fr^R e " t ' d cos £ ip b = 
Tp'b + ih (v ■ fe — w) ip' b — +E b tp b + hQ. R e I0d cos C, ip a = , 



where the prime refers to a derivation with regard to 
the phase £ and v = p/m denotes the center-of-mass 
velocity of the atom with the non-relativistic energy 
£ = p 2 /(2m). Similar techniques have been applied 
more recently also to the Mott scattering of an electron in 
the presence of intense single-mode laser fields [44( , both 
within the relativistic and non-relativistic regime. For 
non-relativistic neutrons, moreover, such an approach 
was taken in Ref. [45| . 

Owing to the (large) mass m of the atom, whose rest 
energy mc 2 is much larger than the photon energy hut, we 
typically have h 2 k 2 / (2mhu>) < 10~ 10 even for ultraviolet 
frequencies and can hence make use of the EA (4(|. In 
this approximation, we will drop the first terms in the 
last system of equations and rewrite them in the form 



ij)' a + iaip a — ifle ?< ^ d cos^ipb, 
4i + %^ h =iQ,e i ^ cosC^a, 



by introducing the following notations 



h (vk — uj 



,13 



n 



h (vk — uS) ' vk uj 



(Al) 
(A2) 

.(A3) 
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The denominator v ■ k lu hereby illustrates the Dopplcr 
shifted radiation frequency as seen by the moving atom 

Eqs. (|A1|) and (|A2|) describe the evolution of the prob- 
ability amplitudes in EA. To find solutions for these two 
first-order equations, we may use the ansatz 



tpb 



to bring them into the simpler form 

A' = ine-^ d cosCe^-W B, 
B' = ifie 1 ^ cosCe- l(Q - /3K ^, 



(A4) 
(A5) 



(A6) 
(A7) 



in which the second terms on the left-hand side of Eqs. 
(|A1[) and (|A2[) have been eliminated. As we will see be- 
low, this re-definition (|A4[) and (|A5P of the probability 
amplitudes facilitates the integration of Eqs. (|A6|) - (|A7I) . 
In addition, we could also decouple these two equations 
by taking the second derivative with regard to the phase 

c, 

A" + (-i(a- /3) + tan C) A' + fl 2 cos 2 C A = (A8) 
B" + (i(a-j3) + tan C) B' + Q 2 cos 2 C B = , (A9) 



with 

Mi,2 = l -{a-P + l±8) , (A14) 

S = ^/(a-p + lf + W , (A15) 

and where A (0) and B (0) refer to "initial" conditions 
with regard to the laser phase, £ = 0. These initial condi- 
tions are fulfilled, for instance, at the origin r = 0, t = 
of the space and time coordinates. Owing to the existence 
of the phase (, however, the solutions (|A12[) and (|A13[) 
are more general and can be analyzed in order to explore 
the time- and space-dependency of the probability am- 
plitudes explicitly, if one wishes to take into account the 
wave vector of the laser field and the motion of the atom 
as a whole. For the atom at rest (p = 0) and if we also 
ignore the fc-dependence of the laser beam, i.e. simply 
perform a replacement C —u>t, the solutions (|A12[) and 
(|A13|) simplify to 



A(t) = \A(0) 



A 2 t Ai . hot 

cos i- — sin 

2 A 2 2 



A 2 t 



i^e- i4 *B(0)sm^ > 
A 2 2 



and for which solutions are known in terms of hyperge- 
ometric functions [22| . However, in this work we are in- 
terested in laser frequencies which are nearly resonant to 
atomic transition frequency. Therefore, we apply RWA 
to obtain solutions to Eqs. (|A8|) and (|A9[) in terms of 
elementary functions. 

In our further discussion, as usual, we consider frequen- 
cies (of the radiation field) which are in resonance with or 
nearly resonant to the atomic excitation, hui ~ E a — Eb- 
In such a resonance regime, it is justified to apply the 
RWA for which the exact solutions can be found for 
Eqs. (|A6|) and (| A7|) . If the counter- rotating terms pro- 
portional to cxp [±i (a — /3 — 1) £] arc ignored on the 
right-hand side, these equations take the form 

A , = M-i<f, de i( a -0+i)C B; (A10) 

B' = —e i +'e- i ( a - f,+1 KA. (All) 
An (exact) solution for A and B is then given by 



( M 1-M2M(C) = - 



^4(0)M2 - ye-***B(0) 



in 



( / i 1 - M2 )s(C) = 



M0)»i - ye-^i?(0) 



B(0)Hi + ye*A(0) 



(A12) 



B(0)fi 2 + y e^M(0) 



, (A13) 



B(t) = tB(0) 



A 2 t .Ai . A 2 t 

cos h 1-. — sin 

2 A 2 2 



Mr %4>d , . . A 2 t 'l _ 
i-r— e Vd B(0) sin -— )■ e 
A 2 2 



huj and A\ 



A\ 



Q 2 R , and are 



with UA\ = E a - E b 
in full ag reement with standard texts on the two-level 
atom [23|]. Using the Eqs. (|A12| and (|A13j) . moreover, 
the conservation of the overall probability of the atom, 
namely of being in one of the two states, |i/> e | 2 -|-|i/> s | 2 = 1, 
can be verified quite easily. 

Although an oscillation of the probability amplitudes 
ip a and tpb can be seen already from ansatz (|A4[) - (|A5[) . 
further insights arc obtained if we specify the "initial" 
conditions for these probability amplitudes and modulate 
the field in a resonance regime with the atomic transition 
frequencies. Therefore, two remarks are in order here be- 
fore we shall further examine the solutions (|A12p - (|A13|) 
for the amplitudes A(Q and -B(C)- First, the initial con- 
dition for C — should be chosen properly. If we assume 
the atom initially to be in the upper state, we have 



A(0) = 1, B(0) = 0, 



(A16) 



for the interaction of the atom with linearly polarized 
field [23||. Second, we still have some freedom in general 
of how to choose the physical parameters, such as the 
momentum p of the atom, its energies E a , Eb of the up- 
per and lower states, the frequency uj of the coupling field 
as well as the Rabi frequency flu. Apart from a suitable 
choice of the two-level atom, these parameters are often 
controlled by the intensity and the propagation direction 
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of the external field(s) acting upon the atom. For ex- 
ample, we may readily fulfill the resonance condition by 
assuming a field with frequency ui, so that 

(1) • Th \ 
1 H J , 

and where n is the unit vector along the propagation 
direction of the field as used before. This resonance con- 
dition gives rise to the simple relation 

a- /3 + l = i) (A17) 

for the reduced quantities (|A3D . 

If we substitute the initial condition (|A16[) into the 
solutions (|A12I) - (|A13|) for a linearly polarized field, the 
p/iase-dependent probability amplitudes read as 

5 5 ' 

B = -\ e^e"^ - 4 e#d e_AtlC > 
2<5 28 

and, together with the resonance condition (|A17|) . take 



the simple form 

A = cos ^ , (A18) 
B = ie**" sin — , (A19) 

where we have used the relations (|A14|) - (|A15|) . Hence, if 
we use the ansatz (fl?) and Eqs. (|A4[) - (|A5|) and then re- 
cover the corresponding expressions of generic notations 
Qr, <j>d and Cj we wm arrive to the form (fT5|) and (|2ll|) for 
the laser-driven two-level atom in collinear- and crossed- 
beam scenario, respectively. 

To summarize this derivation, we have found exact 
solutions (|A18|) - (|A19|) for a particular coupling of the 
atomic motion to the radiation field by applying the 
cikonal, long- wave and rotating- wave approximations. In 
practice, these three approximations arc well justified if 
the frequency of the electric field is nearly resonant to 
the two- level excitation energy (for the RWA), the laser 
wavelength is larger than the atomic sizes (for the LWA) 
and since the de' Broglie wavelength of an atom is typi- 
cally much smaller than the wavelength of the radiation 
field (for the EA). Of course, the RWA might fail if the 
resonance condition is not fulfilled to a sufficient degree 
which may lead to the standard Bloch shift [47| or some 
generalized Bloch-Sicgcrt shift (48j . 
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